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Abstract 

The large availability and rich spectral coverage of today’s observational data 
of the solar corona, and the high spatial and temporal resolution of many 
instruments, has enabled the evolution of three-dimensional (3D) physical 
models to a great level of detail. However, the 3D information provided by 
the data is rather limited as every instrument observes from a single angle 
of vision, or two at the most in the case of the STEREO mission. Two 
powerful available observational techniques to infer detailed 3D information 
of the solar corona from empirical data are stereoscopy and tomography. In 
particular, the technique known as differential emission measure tomography 
(DEMT) allows determination of the 3D distribution of the coronal electron 
density and temperature in the inner corona. This paper summarizes the 
main technical aspects of DEMT, reviews all published work based on it, 
and comments its future development and applications. 
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1. Introduction 

The solar corona can be observed in the white light, EUV, X, and ra¬ 
dio wavelengths. Being the corona optically thin in these spectral ranges, 
its images are two-dimensional (2D) projections of the 3D emitting struc¬ 
ture. Detailed knowledge of the 3D distribution of the fundamental plasma 
parameters of the solar corona (B, N^, Tf) is highly desirable to advance 
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its modeling. Stereoscopy and tomography are powerful observational tech¬ 
niques of the corona, allowing to infer quantitative 3D information of it. An 
excellent introduction to solar stereoscopy can be found in Inhester (2006). 
A recent general review on both techniques covering all the spectral ranges 
listed above can be found in Aschwanden (2011), with a strong focus on stere¬ 
oscopy. In this review we specihcally focus on differential emission measure 
tomography (DEMT) in a more extensive fashion, updating on all published 
work in the held at the moment of writing this article. 

Minnaert (1930) originally developed the scattering theory of the pho- 
tospheric white light (WL) by the free electrons of the corona, that allows 
to infer the 3D distribution of the electron density of the corona from WL 
images, van de Hulst (1950) were the hrst to perform a global corona recon¬ 
struction using eclipse images and assuming full azimuthal axi-symmetry, 
an assumption hrstly relaxed by Leblanc et al. (1970). It was Altschuler 
and Perry (1972) who developed the hrst actual solar rotational tomography 
(SRT) using coronagraph data. A good review on WL SRT can be found in 
Frazin (2000) and Frazin and Janzen (2002), who developed a robust, regular¬ 
ized, positive method for tomographic inversion of the coronal density from 
time series of WL images. Later on, Frazin et al. (2005) hrst introduced 
the concept of DFMT, a technique which uses time series of FUV images to 
determine the 3D distribution of the coronal local-DEM (or LDEM). 

DEMT was developed by Frazin et al. (2009), and hrstly applied by 
Vasquez et al. (2009) to study the 3D structure of coronal prominence cav¬ 
ities. The technique consists of two steps. In a hrst step the tomographic 
inversion of time series of full-sun EUV images is performed, to hnd the 3D 
distribution of the EUV emissivity in each hlter band of the telescope. In a 
second step the emissivities found for all bands in any given coronal location 
are used as a constraint to infer the LDEM. Finally, moments of the LDEM 
are taken, as a result of which global 3D maps of the coronal electron density 
and temperature are produced. 

In this review we summarize the main aspects and applications of the 
DEMT technique. Sections 2 and 3 describe and illustrate the two steps of 
DEMT, section 4 is a review of published results using the technique, and 
section 5 summarizes its main characteristics and future prospects for its 
development and application. 
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2. The Tomographic Model of the Corona 

To perform the EUV tomography, the inner corona volume in the height 
range 1.00-1.25 Rq is discretized on a 25x90x180 (radial x latitudinal x 
longitudinal) spherical grid. Due to optical depth issues (analyzed in detail 
in Frazin et ah (2009)) and EUV signal-to noise levels (which depend on the 
particular hlter), the results are reliable typically in the height range from 
1.03 to 1.20 Rq. 

For each filter band of the EUV telescope separately, time series of full- 
sun EUV coronal images covering a complete solar rotation are used to find 
the 3D distribution of an emissivity-type quantity known as the filter band 
emissivity (FEE). The FEE of each EUV filter is the integral over wavelength 
of the coronal spectral emissivity multiplied by the passband of the hlter. The 
intensity in each pixel is a line-of-sight integral of the FEE. The intensities 
of all pixels of all images can be arranged in a single very large column 
vector, as well as the FEE in every cell (or voxel) of the tomographic grid. 
In this way, both vectors are linearly related through a very large non-square 
sparse projection matrix, that depends on the geometry of the observations. 
Eoth the projection matrix and the pixel intensity vector are known, and 
the problem is to hnd the FEE vector. This poses a non-invertible linear 
problem for each band separately, which is the tomographic problem. 

In the case of both the instrument Extreme ultraviolet Imaging Telescope 
(EIT) on board the Solar and Heliospheric Observatory (SoHO), and the 
instrument Extreme Ultra Violet Imager (EUVI) on board the Solar TEr- 
restrial RElations Observatory (STEREO), the number of EUV bands that 
can be used for DEMT is 3. In the case of the Atmospheric Imaging Assem¬ 
bly (AIA) on board the Solar Dynamics Observatory (SDO) the number of 
bands is increased to 6. 

The 3D distribution of the FEE is found by solving a global optimization 
problem, and the FEE distribution that best reproduces all intensities is 
determined. A thorough technical explanation of all aspects of the inversion 
can be found in Frazin et al. (2009), and discussions on the uncertainties 
involved can be found in Vasquez et al. (2009, 2010, 2011). 

For the EUVI instrument. Figure 2 shows a summary of the EUV tomog¬ 
raphy step. The hrst column shows (from top to bottom) EUVI images in the 
171, 195 and 284 A bands. These images are just one sample from the time 
series actually used to perform the tomography. For each band, the second 
through fourth columns show projected spherical cuts of the tomographic 
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FBE at 1.035, 1.085, and 1.135 Rq, respectively. The last column shows 
the respective synthetic images calculated by integrating the tomographic 
models along the line-of-sight. The black streaks seen in the reconstructions 
near some of the active regions are artifacts caused by the Sun’s temporal 
variability. 



Figure 1: A summary of the EUV tomography. First column shows (from top to bottom) 
EUVI-A images in the 171, 195 and 284 A bands taken near 08:00 on 28 April 2008. 
For each band, the second through fourth columns show projected spherical cuts of the 
tomographic FBE at 1.035, 1.085, and 1.135 Rq, respectively. The last column shows the 
respective synthetic images calculated by integrating the tomographic models along the 
line-of-sight. The black streaks seen in the reconstructions near some of the active regions 
are artifacts caused by the Sun’s temporal variability. From Vasquez et al. (2009). 

To evaluate the accuracy of the tomographic model, the synthetic images 
can be quantitatively compared to the corresponding data images. An ex¬ 
ample of this is shown in Figure 2 from a tomographic reconstruction of the 
solar corona for the bands 171, 193, 211, and 335 A (from tom to bottom) 
of the AIA telescope. 

For each pair of images, the relative difference between the synthetic and 
observed values is below 0.1, 0.2 and 0.3 for 34, 59 and 75% of the pixels. 
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Figure 2: Comparison of synthetic data derived from the tomographic model against 
observed data for the bands 171, 193, 211, and 335 A (from tom to bottom) of the AlA 
telescope. Left panels: data images and corresponding synthetic images computed by LOS- 
integration of the tomographic model. Right panel: frequency histogram of the synthetic 
to observed intensity ratio for every corresponding pair of pixels in the two images. 
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respectively, ±4% depending on the band. The same level of agreement 
holds for off-limb or on-disk pixels considered separately. The tomographic 
model provides then a quite detailed reliable description of the average global 
corona during the reconstructed period. The black rings in the images shown 
in Figure 2 correspond to pixels with projected radius in the range 0.98 to 
1.025 Rq. This near-limb data is not actually used for the tomographic 
inversion, as the emission along their corresponding line-of-sights can affected 
by optically thick emission (Frazin et ah, 2009). 

EUV tomography can currently only be applied from only one or two (in 
the STEREO era) point-of-view. With such limited simultaneous information 
the temporal resolution of the technique is of the order of half solar rotation 
(or about two weeks). Of course, this is the most important limitation of 
the technique, which is then suitable for studying structures that are stable 
during their observed transit. 

3. The 3D Distribution of the DEM 

Once the tomographic step is completed, the FEE of all EUV bands is 
known at each tomographic computational cell. Within each tomographic 
voxel the plasma is expected to be multi-thermal. The LDEM is a measure 
of the thermal distribution within the voxel. 

As each FEE represents emission at a different temperature, they provide 
constraints on the LDEM. Using the passband function of each band of the 
EUV telescopes, and the atomic database CHIANTI version 7.1 to model of 
coronal emissivity. Figure 3 shows the temperature response function (TRF) 
of all EUVI and 4 AIA coronal bands. The temperature range modeled by 
the LDEM is determined by the temperature range of sensitivity of the bands 
that are used. This can be estimated from the FWHM of the TRFSs around 
their respective main peaks. For EUVI and the 3 lower temperature bands 
of AIA, the temperature range is from about 0.5 und up to about 3.0 MK. 
In the case of AIA, this range is expanded up to ~ 4 MK when adding the 
335 A band (see Nuevo et ah (2015)). 

Due to the limited number of available data points (FEE values), and also 
to the narrow band nature of the TRF of each band, the inversion of the 
LDEM function is under determined and implemented by modeling it with 
a family of functions depending on a vector of a small number of parameters 
(typically 3 to 5). When DEMT is based on 3 bands, such as provided by the 
EUVI telescope, the LDEM is modeled by a single normal function. With the 
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Figure 3: Temperature response function of all EUVI coronal bands (dashed, 171, 195, 
284 A), and 4 AIA coronal bands (solid, 171, 193, 211, 335 A). From Nuevo et al. (2015). 

AIA telescope, more coronal bands can be nsed (np to 6), and the LDEM can 
be modeled with combinations of normal fnnctions. The reader is referred to 
Nuevo et al. (2015) for a detailed study on parametric models of the LDEM 
using both telescopes. 

In each tomographic cell independently, the problem consists of Ending 
the values of the parameter vector to best predict with the LDEM the to¬ 
mographic values of all FBEs in that cell. To do so an objective function 
is defined, that measures the quadratic differences between tomographically 
determined FBEs and the those synthesized from the modeled LDEM. 

Once the LDEM is determined, the average electron density in the to¬ 
mographic cell can be computed from its zeroth moment, specifically as the 
square root of the integral of the LDEM over temperature. As an example, 
the left panel in Figure 4 shows a latitude-longitude map of a spherical shell 
of the electron density obtained at one sample height of the tomographic grid, 
for Carrington rotation (CR-)2069, a period of very low magnetic activity. 
At 4 tomographic voxels indicated as I, II, III, and IV in the left panel, the 
right panel shows the normal LDEM model that is found. The error bars 
represent the uncertainty due to the regularization level of the tomographic 
reconstructions. The first moment of the LDEM, divided by its zeroth mo- 
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Figure 4: Left: Latitude-Longitude map of the tomographic electron density Ne in units of 
10®cm“^, at a height 1.075 R©, for the period CR-2069. Right: LDEM at the tomographic 
voxels I, II, III and IV indicated in the left panel. Reproduced from Frazin et al. (2009). 

ment, allows computation of the mean electron temperature predicted by 
DEMT. Examples of spherical shells of the electron density and mean tem¬ 
perature at two different heights of another tomographic reconstruction are 
reproduced in Figure 5. Contour levels of the magnetic-strength 5 of a po¬ 
tential field source surface (PFSS) model are over-plotted, along with the 
magnetically open/closed boundary (see caption in Figure 5). 

4. Results 

Applications of DEMT have included the observational study of coronal 
structures, the use of tomographic results to validate coronal models, and its 
combination with coronal extrapolations of the photospheric magnetic held. 
Following we summarize all peer-reviewed published work on DEMT. 

Vasquez et al. (2009) produced the hrst observational 3D analysis of stable 
coronal prominence cavities, measuring the density and temperature contrast 
between the plasmas in the cavity and in the surrounding helmet streamer. 
As it is characteristic of tomography, their study did not require any ad-hoc 
modeling, as needed in the forward modeling approach. 

Being suited to study coronal structures that are stable over half solar 
rotational time, tomography works best at solar minimum. Vasquez et al. 
(2010, 2011) analyzed the global thermodynamical of the solar corona during 
the minimum of activity between solar cycles 23 and 24, and discussed their 
results in relation to the fast and slow components of the solar wind. These 
works also include comparisons to other observational non-tomographic stud¬ 
ies of the same periods, providing cross validation results for the technique. 
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Figure 5: An example of the 3D reconstruction of the thermodynamical state of the solar 
minimum corona, during the period CR-2077. Latitude-longitude maps of the electron 
density (top) and mean electron temperature (bottom) at heights 1.035 (left) and 1.075 
Rq (right), as derived with the DEMT technique. Solid-thin curves are magnetic strength 
B contour levels of a PFSS model of the coronal magnetic field, with white (black) rep¬ 
resenting outward (inward) oriented magnetic field. The solid-thick black curves indicate 
the location of the magnetically open/closed boundary. 


As an example of a 3D reconstruction of the solar minimum, Figure 5 dis¬ 
plays latitude-longitude maps of the electron density and mean temperature, 
derived with the DEMT technique for the period CR-2077. At both heights, 
the location of the open/closed boundary of the PFSS model is characterized 
by a very high transverse gradient in both the electron density and the mean 
temperature maps derived from the DEMT analysis. 

DEMT can be combined with the global PFSS magnetic models of the 
solar corona, an approach dubbed the Michigan Loop Diagnostic Technique 
(MLDT). DEMT results are traced along the held lines of the magnetic 
model, allowing study of the thermodynamical properties of magnetic hux 
tubes in the quiet sun in a statistical fashion. Huang et ah (2012) applied 
MLDT to study one rotation during the last solar minimum. They found 
the ubiquitous presence of magnetic loops with downward gradients of tem¬ 
perature dubbed down loops. Down loops were found to be dominant in the 
latitude range ±30° (left panel in Fig. 6). Nuevo et ah (2013a) extended 
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Figure 6: Left: Location of loops with negative (down) and positive (up) temperature 
gradient with height, in blue and red, respectively, for CR-2077. Right: Evolution of 
the fraction of down (up) loops indicated as light-blue (red) diamonds, and the sunspot 
monthly number divided by 32.1 (dark-blue squares). The correlation between the fraction 
of down loops and the sunspot number is indicated. Adapted from Nuevo et al. (2013a). 


the study to a sequence of rotations that included the solar minimum. Their 
study revealed a clear anti-correlation between the global coronal activity 
level and the number of down loops present in the corona (right panel in 
Fig. 6). They found that down and up loops are characterized /3 ~ 1 and 
/S < 1, respectively, and proposed an interpretation of their results in terms 
of Alfven wave damping. 

DEMT results have been also used to constrain and validate MHD models. 
Reconstructions of the electron density and temperature have been used as 
a constraint to 3D MHD solar wind models coupled to the Space Weather 
Modeling Framework (SWMF) code suite, that solves for the different elec¬ 
tron and proton temperatures (van der Holst et ah, 2010). Improvements 
on the performance of the models include a more accurate prediction of the 
occurrence and density of co-rotating interaction regions in the heliosphere. 

Global two-temperature models of the corona and inner heliosphere have 
been validated in the inner corona with DEMT reconstructions of the elec¬ 
tronic plasma parameters (Jin et ah, 2012), with most of the model outputs 
fitting the observations very well, as seen in Figure 7. The 3D products of 
DEMT have also been used as a validation tool for a study on coronal heat¬ 
ing by surface Alfven wave damping, implemented within the MHD model 
of the solar wind in the SWMF (Evans et ah, 2012). DEMT results have 
been recently used as a validation tool in a study of the charge state com¬ 
position of the slow solar wind derived from an ionization evolution code 
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Figure 7; DEMT validation of the electron density (left) and temperature (right) of a 
two-temperature MHD model of the solar corona, for the period CR-2077. The color 
rings spanning the height range 1.0 to 1.25 Rq show, in a sample meridional cut of the 
corona, the MHD to DEMT ratio of the respective quantity. The outer greyscale maps 
(r > 1.25 Rq) displays the corresponding MHD results. Reproduced from Jin et al. (2012). 


coupled to the wave-driven MHD solar wind model of the SWMF (Oran 
et ah, 2014). DEMT allowed detailed validation of the latitudinal transition 
of the electronic plasma parameters in the open field lines surrounding the 
coronal equatorial streamer belt (Figure 8). 

All DEMT studies reviewed so far were based on EUVI data, having 3 
EUV bands with a sensitivity range ~ 0.60 — 2.70 MK. More recently, Nuevo 
et al. (2013b, 2015) extended the DEMT technique to use the 4 cooler AIA 
bands (aimed at studying the quiet sun), sensitive to the range ~ 0.55 — 3.75 
MK. Their study corresponds to CR-2099, a rotation of the rising phase of 
the current solar cycle 24. While in previous studies the LDEM was always 
modeled by a single normal distribution, the extra AIA hlter, and the in¬ 
creased sensitivity range, allowed exploration of new parametric LDEM mod¬ 
els. Using 4 bands, the model that consistently achieves the best predicted- 
to-reconstructed FBEs is a bimodal distribution, being a superposition of 
two normal distributions with distinct cool and hot components. The mean 
centroids of the two components in the quiet diffuse corona are found to be 
logio (^o,i) = 6.15 and logio {To, 2 ) = 6.42, values that are very consistent 
with independent determinations of the characteristic temperatures of the 
solar corona, as discussed in the same study. 
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Figure 8: DEMT validation of a 3D wave-driven MHD model of the solar wind. Electron 
temperature versus angular distance from the magnetically open/closed boundary (aver¬ 
aged over all longitudes) for the period CR-2063, derived from the MHD model (red) and 
from the DEMT inversion (black). Reproduced from Oran et al. (2014). 


The square electron densities of the two components are and N ‘^21 
respectively, so that the total square electron density of the LDEM is iVg = 
Afg 1 + A measure of the bi-modality is then the ratio (Afe, 2 /Ae)^- The 
study of this ratio throughout the diffuse quiet corona reveals that the bi- 
modality of the LDEM is ubiquitous, and that it is stronger for denser regions, 
as shown by the left panel of Figure 9. 

Nuevo et al. (2015) also validate the LDEM inversion technique by ap¬ 
plying it to standard 2D DEM studies. Examples of the bimodal DEMs are 
shown in the right panel of Figure 9. The DEMT study shows that LDEM 
of the quiet corona is bimodal at the spatial resolution of the tomographic 
grid, which is 0.01 Rq x 2° x 2°, or about (7 x lO^km) x [(2.44 x 10"^km)^] 
for a representative voxel at a height of 0.1 Rq above the photosphere at 
the equator. The authors argue that the nanoflare heating scenario is less 
likely to explain these results, and that alternative mechanisms, such as wave 
dissipation appear better supported by them. 
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Figure 9: Left: Scatter plot of the measure of the bi-modality (A^e, 2 /.^e)^ (see text) versus 
Nf. derived with AIA, at 1.075 Rq. Right: Bi-modal DEM curves obtained for selected 
regions within and around an active region. Red/yellow curves correspond to hotter and 
denser regions, blue/dark curves correspond to colder and less dense regions. Reproduced 
from Nuevo et al. (2015). 

5. Concluding Remarks and Future Prospects 

DEMT provides a quantitative average description of the solar corona over 
a full solar rotation with the following main characteristics, 

• The products are 3D maps of: a) the hlter band emissivity in each 
EUV band, b) the local-DEM, and c) its moments, such as (N^) and 
(Te). 

• The temporal resolution of the technique is currently limited to the 
transit time of coronal features, i.e. ~ 14 days. DEMT reconstructions 
are then reliable descriptions of slowly evolving coronal structures, such 
as the diffuse quiet corona, and coronal holes. DEMT is not suited to 
study fastly evolving ARs. 

• The spatial resolution depends on the computational grid size, which 
in turn is constrained by the cadence of the image series. Typically, one 
image every 6 hrs is used, a time over which the Sun rotates about 3.3°. 
The adopted tomographic cell size is then 2° in both angular directions 
and 10“^ Rq in the radial direction. 

• The image processing, tomographic inversion, and DEM determination, 
are fully automated tasks, with little user interaction. 


13 






• Its implementation does not require any ad-hoc modeling. 

• Finite FOV and finite computational grid effects, issues of relevance for 
white light tomography, are not important for EUV tomography. 

The current implementation of the DEMT and MLDT techniques involves 
a suite of codes written in the C and IDL programming languages, which 
can be efficiently ran in nowadays desktop multi-core computers. Future 
immediate planned applications of DEMT and MLDT involve comparative 
studies of the last two solar minima (using the EIT and EUVI instruments), 
as well as the study of the coronal radiative losses predicted by MLDT in the 
quiet sun. 

Tomographic inversion of time series of full-sun images is the only available 
observational technique that can provide global constraint and validation 
to large scale MHD modeling of the corona and the solar wind. As such, 
DEMT is a highly valuable tool to help in the continued development of 
global coronal models. DEMT could provide highly valuable 3D maps of the 
detailed temperature distribution of the coronal plasma if full-sun spectral 
images (of the type provided by the EIS instrument) become available. It 
would be highly desirable that such an instrument will be operational during 
the a solar minimum period, when single point-of-view tomography is most 
fruitful. 

The authors thanks CONICET grant PIP lU Nro 11420100100151 to 
lAFE that has funded this research. The author also thanks the collabo¬ 
ration and useful comments by Federico A. Nuevo and Richard A. Frazin. 
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